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■^^ ' Abstract 

OO ' Let r be a random field invariant under the action of a compact group G. In the line of previous 

work we investigate properties of the Fourier coefficients as orthogonality and Gaussianity. In 
particular we give conditions ensuring that independence of the random Fourier coefficients 
implies Gaussianity. As a consequence, in general, it is not possible to simulate a non-Gaussian 
invariant random field through its Fourier expansion using independent coefficients. 
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^ : 1 Introduction 

in 



Recently much interest has been attracted to the investigation of properties of random fields on the 
sphere §^ that are invariant (in distribution) vi^ith respect to the action of the rotation group 50(3), 
highlighting a certain number of interesting features (see [9], [1] e.g.). This interest is motivated 
^^ . mainly by the modeling and the investigation of cosmological data. 

For instance in [1] it was proved that assumptions of independence of the Fourier coefficients of 
the development in spherical harmonics 



T = ^ ^ atmYlm (1-1) 



of the real random field T , in addition to invariance, imply Gaussianity. More precisely it was proved 
that if the coefficients Ohn, ^=l,2,...,0<77i<^ are independent, then the field is necessarily 
Gaussian (the other coefficients are constrained by the condition a^,-m — (— l)™a^„). This result 
implies, in particular, the relevant consequence that a non Gaussian invariant random field on S^ 
cannot be simulated using independent coefficients. 

It is then a natural question whether this property also holds for invariant random fields on 
more general structures. A result in this direction was obtained in [2] where it was proved that 
a similar statement holds in general for an invariant random field on the homogeneous space of a 
compact group, provided the development is made with respect to a suitable Fourier basis satisfying 
a particular condition (see Assumption 4.3 below). 



The object of this paper is to pursue the hne of investigation started in [2], in particular in the 
direction of determining in which situations Assumption 4.3 is true. Actually we do not know 
whether this condition holds in general and its investigation is one of the main objects of this paper. 
Until now it is known to hold in the case of §^ for the basis given by the spherical harmonics (Yfm)m, 
iU>l (see [9] p. 145). 

We shall give a new condition, equivalent to Assumption 4.3, which will enable us to prove that 
it is satisfied for every self-conjugated Fourier basis of the sphere §^ (and not just the spherical 
harmonics) and also for an important class of self-conjugated bases of the sphere S'^. In particular 
it is non possible to simulate a non Gaussian invariant real random field on §^ using independent 
coefficients with respect to any self-conjugated bases of the irreducible G-modulc of L^(S^). 

Besides this characterization of Gaussianity, we discuss other properties of the Fourier coefficients 
of an invariant random field as orthogonality (it is well known that the ae^s of the development 
(1.1) on §^ for an invariant random field with finite variance are known to be orthogonal, see [9], 
p. 140) and invariance of their distribution with respect to rotations of the complex plane. 

The plan of the paper is as follows. In §2 we recall the main properties of the Fourier development 
of a random field on the homogeneous space <!?^ of a compact group and give necessary and sufficient 
conditions for its invariance in terms of its development. In §3 we investigate properties of its 
coefficients as orthogonality, among other things. It turns out that, unlike the case of S^, they 
are not orthogonal in general and precisions are made concerning this phenomenon. In §4 we recall 
(from [2]) results giving the characterization of Gaussianity which is our main concern. As mentioned 
above, these results, in many cases of interest, hold under the assumption that the Fourier basis that 
is chosen for the Fourier development enjoys a certain property (Assumption 4.3) with respect to 
the action of the group. This section also contains a converse result, giving conditions on Gaussian 
coefficients in order to produce an invariant random field. 

The remainder of the paper is devoted to the investigation of the validity of Assumption 4.3. 
In §5 we give a new equivalent condition which is the main tool for the investigation of the two 
main examples (§^ and S'^) which are the objects of §6 and §7. We are also able to prove that 
independence of the Fourier coefficients entails Gaussianity in some situations of interest in which 
it is known that Assumption 4.3 does not hold (Theorem 6.4), in particular covering the case of the 
basis of the spherical harmonics {Yirn)m for i? = 1 (which completes the proof of the result of [1]). 

Finally §8 points out some open questions. 

2 A.s. square integrable random fields and Fourier develop- 
ments 

Throughout this paper ^ — GjK denotes the homogeneous space of a compact group G, K being 
a closed subgroup. We denote x — t- gx, g E G the action of G and dg and dx respectively the Haar 
measure of G and the G-invariant measure on .^ . We assume that both dx and dg are normalized 
and have total mass equal to 1 and we write L'^{^) for L'^{^,dx) and L'^{G) for L^{G,dg). The 
spaces L^ are spaces of complex valued square integrable functions. 

Let us denote by G the set of equivalence classes of irreducible representations of G and let, for 
every a G G, H" a Hilbert unitary G-module of class a fixed from now on. For every / G L^{G) let 

7(a):=Vdh^ //(g)a(g-i)dg. 

JG 

f{a) is a linear endomorphism of H"^ . The Peter- Weyl theorem (see [12] or [3] e.g.) states that 

fig) =. J2 ^/dh^ tr(7(a)a(g)) . (2.1) 

aed 



It is immediate that this is a good definition not depending, in particular, on the G-invariant scalar 
product that is considered on the G- module W^ . 

As soon as an orthonormal basis hi, . . . ,hm, in = dima, of if^ is chosen, one can define the 
matrix coefficients 

and the corresponding matrix entries of /(cr): 

f{a)^j = Vdh^ f f{g)D^g-^)dg = Vdh^ / f{g)DfJ^dg = Vdh^(/,i^J,)i2(G) 
Jg jg 

so that (2.1) becomes 



dim (T dim <j 

fig) ^ J2 ^/di^ E J^'^'^h DUg) - Y. d™^ E (/'^J.)i^(G) D^g) (2.2) 

CTeG *J = i aeG «J = i 

The Peter- Weyl theorem also entails that, for every orthonormal basis hi,... ,h„i, rn — dinicr of 
H'^, the normalized matrix elements Vdim a Df, , cr e G, 1 < i,j < dima, form an orthonormal 
complete basis of L'^{G), so that (2.2) is just the corresponding Fourier development and ,f{a)ij is 
the coefficient corresponding to the element Vdim a D'J^ of this basis. 

Let us denote Lg the left action of G on L'^{G), that is Lgf{h) — f{g^^h). It is immediate that 
the functions (-D'^(.g)ij)i<i<dim<T, appearing in the columns of D"^ , span a subspace of L'^{G) that is 
invariant and irreducible with respect to this left action. 

From (2.1) a similar development can be derived for i^(^), ,9^ = G/ K. Actually remark that, 
for a fixed xq G ^, the relation f{g) ~ f{gxa) uniquely identifies functions in i^(^) as functions 
in L^ (G) that are right invariant under the action of the isotropy group K of xq (that is formed by 
the elements k € G such that kxo = xq). For such functions / we have, for every k € K, 



f{<^) = / I{g)<y{9-')dg = / f{gk)a{g-')dg - / f{t)<y{kt-')dt = a{k) / f{t)a{t-')dt = a{k)f{a) 

JG JG JG JG 

which implies that /(ct) is -ff^-valued, H^^ denoting the subspace of H'^ of vectors that are invariant 
under the action of K, i.e. /(a) G End{H'^ , H'^) . Hence, for a choice of an orthonormal basis 
hi, . . . , hfn of H"^ such that hi, . . . ,hk span H'^, the matrix f{a) will have all zeros in the rows 
from the (fc + l)-th to the vn-th. For / G LF'{^) we shall consider, for simplicity, that f{(j) is a 
dime X diiuCT matrix with zeros on every row but for the first dimTf^ ones, corresponding to the 
first elements of the basis, that are supposed to be X-invariant. Remark that H"^ might be reduced 
to {0}. 

As a consequence of the aforementioned Peter- Weyl theorem we have the decomposition, that we 
shall need later, 

^'(■^) = y^" (2-3) 

(jgg l<i<dim(_H"g-) 

where G denotes the set of equivalence classes of irreducible representations of G and the V^ are 
irreducible G-modules of class a. 

We consider on .^ a real or complex random field {T{x))x£X ■ This means that there exists a 
probability space (r2, ^, P) on which the r.v.'s T{x) are defined and we shall always assume joint 
measurability, i.e. (a;,cj) -^ T{x,uj) is 3§{,'^) (E) ^ measurable, SS{!^) denoting the Borel a-field of 

jr. 



T is said to be a.s. continuous if the map x — > T{x) is continuous a.s. It is said to be a.s. square 
integrable if 

/ |T(x)pdx < +00, a.s. (2.4) 

Remark that a.s. square integrabihty does not imply existence of moments of the r.v.'s T{x). If T 
is a a.s. square integrable random field on G then the function x -^ T(x,uj) belongs to L'^(^) a.s. 
and one can define "pathwise" 

f(cr) = Vdimcr / T{g)a{g-^)dg (2.5) 

Jg 

which is now a End{H'^ )-vahicd r.v. Similarly we have the analog of (2.1), i.e. 

^(^) ^ Yl ^dimcr tr(f (cr)cr(/i)) (2.6) 

cred 
or 

dim (T dim a 

T{h) = Y. ^di^ E ^(^)*^- ^U^) -^H'^'-^^'ll (^' ^J^)2 ^J.e^) (2-7) 

the series converging a.s. in L'^{G). 

For a random field T we define the rotated random field T^ as T3{x) ~ T{gx). 

Definition 2.1 ^ a.s. square integrable random field T on 1%' is said to he G-invariant if, as a 
L'^{J%') -valued random variable, it has the same distribution as the rotated random field T^ for every 
g G G, in the sense that the joint laws of 

{T{xi),...,T{xm)) and {T{gxi), . . . ,T{gXm)) (2.8) 

coincide for every g £ G and xi, . . . , x„j £ JT. 

More generally a family {Ti)i(^,y of random fields on 2y is said to be invariant if and only if for 
every choice of g (£ G, ii, . . . , im G ^ and xi, . . . , x^ G -^ , the joint laws of 

(rij(xi), . . . ,Ti^(xm)) and (Ti^{gxi), . . . ,Ti^{gxm)) (2.9) 

coincide. 

The following will have some importance later. We thank D.Marinucci and G.Peccati for informing 
us of the existence of this result. 

Proposition 2.2 Let T a a.s. square-integrable invariant random field on 2^ and define, for f e 



T{f):^ / T{x)f{x)dx (2.10) 

Then, for every g £ G and every /i, . . . , /,„ G L'^{G), the two random variables 

(r(/i),...,T(/„)) and (T9(/i),...,T9(/„0) 
have the same distribution. 
Proof (See [10]. 



Proposition 2.3 Let T a a.s. square integrable random field on G. Then T is invariant if and only 
if, for every g Cz G, the two families of r.v's 

are equi- distributed. 

Proof. Let us assume T invariant and let cr G G. Then for every v^w G H'^ the function g — > 
{a{g^^)v^ w) is bounded and therefore in L?'{G). Therefore, according to Proposition 2.2 and denot- 
ing by ^ equaUty in law, we have for every g G G, 

{T{a)v, w) = Vdimcr / T(K)(a(hr^)v, w) dh - Vdima / T{gh){a{h-^)v, w) dh = 
Jg Jg 



= Vdimcr / T{t){a{t-'^g)v,w)dt^\/ dim a / T{t){a{t-^)a{g)v,w) dt ^ {T{a)a{g)v,w) 
Jg Jg 

This being true for every v,w £ H'^ , we have that, as i<^n(i(-ff'^)-valued r.v.'s, T{a) and T{<T)(T{g) 
have the same distribution. Quite similarly, only in a just more complicated way to write, 

(T(cri ),..., r(o-„)) and (T(cri)cri(g), . . . , T(cr„)cr„(g)) 

have the same distribution as a End{H°'^) ® • • • i<^n(i(_ff'^")-valued r.v., thus proving the only if 
part of the statement. The converse follows easily from development (2.6). 

■ 
Let / G L^(^) and V C L^(^) an irreducible G-module. We can then consider the orthogonal 
projection Pyf of / on V. Similarly for a a.s. square integrable random field T on ,^ let us denote 
Ty its orthogonal projection on V. Remark that by definition (the functions of V are necessarily 
continuous) Ty is always a continuous random field. 

Let us denote by Dij {g) the matrix elements of the left regular action of G on F with respect to 
the orthonormal basis {vi)i of V and let us consider the random coefficients of the development of 
T with respect to this basis 



T{x)vi{x)dx. (2.11) 

% 

We denote by a the complex vector with components a^, i — 1, . . . ,(i. Then the coefficients of the 
rotated random field T^ are obtained through the relation 



T{gx)vi{x) dx — / T{x)vi{g ^x)dx = 
s: J a: 

d d (2.12) 



= ^ Dk^{g) f T{x)vk{x) dx^J2 Dk^{g)T{vk) = ^ Afc(ff~' 

fe=l ^ fe=l k=l 

aP = D{g-^)a (2.13) 

d 

Tv{x) = ^akVk{x) , 
fe=i 

it is immediate that Ty is invariant if and only if the random vectors a and D{g)a have the same 
distribution. 



that is 
As 



With respect to the Peter Weyl decomposition (2.3) we have 

dimHJ 

Using the fact that the projectors are G-equivariant (i.e. commute with the action of G) it is easy to 
prove the following, not really unexpected, statement (anyway see Proposition 3 of [11] for a proof). 

Proposition 2.4 T is invariant if and only if the family iTvf)„^Q i<i<dim_w '^■f '''O'^^dom fields is 
invariant. 

We shall therefore concentrate our attention mainly on the projected random fields Ty. When 
dealing with a real random field it is natural to require that the basis wi , . . . , w^ of the G- module 
V , with respect to which the coefficients are computed "respects" the real and imaginary parts and, 
in particular, \iV — V, that this basis is stable under conjugation. As explained in [2], §2 and the 
Appendix, it is actually possible to decompose L^(<^) into a direct sum of G- modules in the form 

L\,T) = y, ® (F, e T^) (2.14) 

where the direct sums are orthogonal and 

i e J^° ^ Fj = Fi, i(:,J+-^ViLVi. 

We can therefore choose an orthonormal basis {vik)ik of L'^{^) such that {di = dim 14) 

• for i G ^° , (vik)i<k<di is an orthonormal basis of Vi stable under conjugation; 

• for i G ^^, {vik)i<k<di is an orthonormal basis of Vi and {u[k)i<k<di is an orthonormal basis 

ofF,. 

It is immediate that if T is a real random field and i G ^° then Ty. is also a real random field. 
On the other hand, if z G J^^ then Ty and Ty may not be real (actually they cannot be real unless 
they vanish) , whereas Ty + T— will be real. 

Remark 2.5 Representations of a compact group G are classically classified as of real, complex or 
quaternionic type (see [3], p. 93 e.g.). In order to be self-contained let us recall that a conjugation 
J of a G-module V is an antilinear {J{av) — aJ{v)) equivariant map J : V — > V. 

A G-module V is said to be real if there exists a conjugation J : V ^^ V such that ,P = \ and 
quaternionic if there exists a conjugation J :V ^f V such that J^ = — 1. It is complex if it is neither 
real nor quaternionic. 

The important thing is that an irreducible G-module is of one and only one of these types and 
that equivalent G- modules are necessarily of the same type. If an irreducible G-module V C L^(^) 
is such that V = V, the usual conjugation J : u — > u is a real conjugation, so that V must be of real 
type. In particular, if a representation is of quaternionic or complex type, it cannot contain in its 
isotypical space a G-module that is self-conjugated, so that in the decomposition (2.14) it cannot 
be of type J^°. 

The irreducible representations of even dimension of SU(2) are quaternionic and the corresponding 
G-modules appearing in the Peter- Weyl decomposition of this group cannot, therefore, be self- 
conjugated. 

3 Properties of the coefficients 

In this section we give results concerning two properties that are enjoyed by the coefficients T{(j)ij , 
c G G, 1 < i,j < dimcr, of the Fourier development of an invariant random field on ^ . 



A random field T is said to have finite variance if 



(/ |r(x)|2rf.T) <+00 (3.1) 



Remark 3.1 If (3.1) holds, then the map x -^ T{x) necessarily belongs to L'^{^) a.s., so that 
T is a.s. square integrable. Also, by the Cauchy-Scwartz inequality, if T has finite variance, the 
random variables T{f), f G i^(<^), defined in (2.10), have finite variance. In particular the Fourier 
coefficients of T, with respect to any Fourier basis, also have finite variance. 

It is well known (see [1], [9] p. 126) that in the case ^ = §^, G — 5*0(3), if T is invariant and has 
finite variance, its Fourier coefficients with respect to the basis formed by the spherical harmonics 
are pairwise orthogonal. Our first concern in this section is to investigate this question in the case 
of a more general basis and for a general homogeneous space of a compact group. 

Remark that, if T is invariant and has finite variance, for every a that is not the trivial represen- 
tation, for the matrix entries T{(j)ij we have 



E[f (a)] = E[f (a)a(.g)] = E[f (a)] / a{g) dg ^ (3.2) 

Jg 



'G 

Theorem 3.2 Let T a finite variance invariant random field on 2y and (Ti, (72 G G. 

a) If (7 1 and 02 are not equivalent, then, for every orthonormal bases of H"^ and H'^^ the r.v.'s 
T{ai)ij and T{ij2)m o^Te orthogonal, 1 < i,j < dinicri, 1 < k,i < dim(T2- 

b) If (Ji ^ 02 ~ <y let T(a) = E[T(cr)T((T)*]. Then Cov {T (a) ij,T (a) ke) = dj£T{a)ik- In particu- 
lar coefficients belonging to different columns are orthogonal and the covariance between entries in 
different rows of a same column does not depend on the column. 

Proof Recah first that T{a)ij = {T,Dj-)2 so that thanks to Remark 3.1 the rv.'s T'(cr)y's have 
themselves finite variance. Let us denote DfAg) the matrix elements of the action of G on H'^ . 

a) As {T{ai),T{a2)) has the same joint distribution as {T{ai)ai{g),T{a2)cr2{g)) for every g & G 
by Proposition 2.3, 



E[Tia,),,Tia2)k,]^E[iTia^)D^I{g))iTia2)Dllig))] = 

dim fTi dim a2 

= E E DrH9)DZi9'Hfi^lhrn<J2)km] . 
r—1 m— 1 

This being true for every (7 e G, it is also true if we take the integral of the right hand-side over G in 
dg. As the functions D^- and D'^^ are orthogonal for every choice of the indices, the representations 
CTi and (72 being not equivalent, we find 



E[T{<Ji)^jT{a2)M] - . 
h) li ai — a2 — (J, the previous computation gives 

dim a 



E[f{a),,f{a)k,] - J2 E[f(a)„f(a)fc„] / D-,^{g)D'^,{g) dg 



& 



im (7 



E E[T(cr)jrT((T)fcm](5rm(5j£ = (5j£ ^J ^\T {u) irT {a) kr] =5jf:'^{(y)ik 



Theorem 3.2 states that the entries of T{a) might not be pairwise orthogonal and this happens when 
the matrix F is not diagonaL This phenomenon is actually already been remarked by other authors 
(see [8] Theorem 2 e.g.). 

Example 3.3 below provides an instance of this phenomenon. Of course there are situations in 
which orthogonality is still guaranteed: when the dimension of H'^ is one at most (i.e. in every 
irreducible G-modulc the dimension of the space H'^ of the ii'-invariant vectors in one at most) as 
is the case for G ~ SO{d), K ~ SO{d — 1), G/K = S''^^. In this case actually the matrix T{a) has 
just one row that does not vanish and T(a) is all zeros, but one entry in the diagonal. 

In order to produce an example of invariant random field for which the entries of r(cr) arc not 
pairwise orthogonal, we recall first a well known definition. 

Let Z — Zi + iZi a complex r.v. Z is said to be Gaussian complex valued if (Zi,Z2) is jointly 
Gaussian. Z is said to be complex Gaussian if, in addition, Zi and Z2 are independent and have 
the same variance. If Z is centered this is equivalent to the requirement that their distribution is 
invariant with respect to rotations of the complex plane. We shall use the following properties. 

• A centered Gaussian complex valued r.v. Z is complex Gaussian if and only if E[Z^] = 0. 

• Two centered complex valued Gaussian r.v.'s Zi, Z2 are independent if and only if 'E[ZiZ2\ — 
E[ZiZ2] =0. 

Example 3.3 Let a E G and V C L'^(G) an irreducible G-module of dimension d > 2 of class cr 
and denote by the matrix D'^{g) the action of G on y with respect to a fixed basis. Let Zi, . . . , Zd 
independent centered complex Gaussian r.v.'s such that E[|Zjp] — 1 for every j. Let B = {bij)ij 
the random matrix defined as bij = ouZj, Ui G C. Then the random field 

is invariant and, as it is immediate that T{a) — B, its coefficients T(a)ij are not pairwise orthogonal. 
Let us check invariance. Let G = T{a)D''{g), then c^ — aiX]fe=i ■^kD'^Ag) = aiWj where 



W,^Y.^^^kji9) 



k=l 

In view of Proposition 2.3 we must therefore just prove that the W^j's are complex Gaussian, inde- 
pendent and that E[|T4^jp] — 1. First it is immediate that they are Gaussian complex valued. We 
have also 

d d 

E[iy,Wfe]=E[^ Zn'z;.Dl^{g)DfJg)\ - ^ 5nrDl^{g)DfJi) = 

h.r—1 h,r—l fo o\ 

d d ^■^■■^' 

r—l r—1 

Similarly, as F,[ZhZr] — for every 1 < h,r < d (recall that E[Z^] == for a centered complex 
Gaussian r.v. Z), 

d 

E[W,Wk] = E [ E Z^,ZrDZ,{g)D'^,{g)] = . (3.4) 

h,r=l 

(3.3) for k = j gives E[|T4^jp] = 1, whereas (3.3) and (3.4) together imply that Wj and Wk, k ^ j, 
are independent. Finally (3.4) for k = j gives E[T4^^] = for every j so that the Wj's are complex 
Gaussian, which completes the proof. 



Corollary 3.4 Let T an invariant random, field with finite variance on .% and let V C L^(^) an 
irreducible G-module different from the constants. Then the coefficients {ak)k of the development of 
the projection Ty of T on V with respect to any orthonormal basis of V are centered, orthogonal, 
and have a common variance c. 

Proof. It is repetition of the arguments of the proof of b) of Theorem 3.2. As pointed out in Remark 
3.1 the coefficients a^'s have themselves finite variance and, thanks to (2.12) and V being different 
from the constants, they are also centered. From (2.13) we have, for every g E G, 

d 
E[afea?] = ^{{D{g)a)},{D{g)a)i\ ^ ^ Dkr{g)Dij{g)E[araj] . 

Integrating in dg and using the orthonormahty properties of the matrix elements Dij (g) we find 

_. dim V .. dim V 

EfafcoZ] = — — — ^ SkiSrjE[araj] = — — 77 fe ^ E[|aj|^] . 
For k ^ i this gives immediately the orthogonality, whereas for k ^ £ wc have 

dimy 

dimy 



E[Kn-;^EE[KI^] 



so that the a^j's have the same variance. 

■ 
Another feature appearing in the case G — 50(3), ^ = S^ is that the coefficients aim of the 
development in spherical harmonics (1.1) of an invariant random field have each a distribution that 
is invariant with respect to rotations of the complex plane ifm^O. The following discussion aims 
to see what can be said in general concerning this property. 

Remark 3.5 Let V C L'^(^) be an irreducible G-module of dimension d > 1 and T C G a maximal 

torus. Let 

d 

V = ^Uk (3.5) 

fc=i 

be a decomposition of V into orthogonal irreducible components of the action of T on V. As T 

is abelian, d[m{Uk) = 1 for every k — l,...,d. Let Uk G f7fc be a unit vector. Then LfUk = 

Ukit^^') = Xk{t)uk for f G T, where Xk denotes the character of the representation of T on Uk- If 

we consider the Fourier development of an invariant random field T with respect to the orthonormal 

basis (ui, . . . ,Ud) 

d 

T = ^ UkUk 
fc=i 

then, as, for t G T, 

d d 

T{t'^x) = ^ akUk{t^^x) = ^ akXk{t)uk{x) 

k=l k=l 

and T(t~^x) and T{x) have the same distribution for every t G T, necessarily, for every k such 
that the action of T on Uk is not trivial (that is Xk{t) ^ 1), the coefficient a^ must be invariant in 
distribution with respect to rotations of the complex plane (and therefore, if it is Gaussian, it must 
be complex Gaussian). 



Remark also that the action of T over V cannot be trivial, that is Xk ^ ^ for some k necessarily. 
Actually, as all maximal tori arc conjugated (that is if T' is another maximal torus then T' ~ gTg~^ 
for some g G G) then the action of all maximal tori on V would be trivial which is impossible as the 
union of all maximal tori is the group itself so that this would imply that the action of G itself is 
trivial, whereas we assumed V to be irreducible and with dimension d > 1. 

The property, mentioned above, of the random coefficients with respect to the basis of the spherical 
harmonics in the case of the §^, appears now as a particular case. 

4 Invariant random fields with independent Fourier coeffi- 
cients 

In this section we see results that state that independence assumptions on the Fourier coefficients 
implies Gaussianity of the coefficients and of the corresponding random field. 

Theorems 4.1 and 4.4 below are already known (see [2]) and we reproduce them only to be self- 
contained, our main concern being the investigation of the validity of Assumption 4.3, which is a 
necessary condition in many situations of interest. 

Theorem 4.1 Assume G to be connected and let T a a.s. square integrable G-invariant random 
field on the homogeneous space J^ . Let V an irreducible G-module of L'^{^) with dimension d > 1 
and let us assume that coefficients {ak)k of the development 

d 

Tv = 2^ °'^'"^ 
fe=i 

with respect to an orthonormal basis {vk)k ofV are independent. Then they are necessarily Gaussian 
and the random field Ty is Gaussian itself. 

Remark that in the statement of Theorem 4.1, as in Theorem 4.4 below, we make no assumption 
concerning the integrability or the existence of finite moments of the r.v.'s T{x) and/or Ufc. But, of 
course, under the assumptions of the theorem it follows that necessarily the r.v.'s Tv{x) and Ok has 
finite moments of every order. 

The proof of Theorem 4.1 relies on the following Skitovich-Darmois theorem, actually proved in 
this version by S. G. Ghurye and I. Olkin [6] (see also [7]). 

Theorem 4.2 Let Ai,...,Aj. be mutually independent random vectors with values in R". If, for 
some real nonsingular n x n matrices Aj, Bj, j — 1, . . . , r, there are two linear statistics 

r r 

ii = 2^ AjXj , L2 = 2^ BjXj 

that are independent, then the vectors Xi , . . . , Xr are Gaussian. 

Proof oi Theorem 4.1. Let us denote again by D{g) the representative matrix of the left action of 
g Cz G on V with respect to the orthonormal basis (ufc)fc and by a the vector of the coefficients a^. 
Thanks to (2.13) we have 

distr „ „. _l^ 

a = a^ — U[g )a . 

Let 1 < fci < ^2 < dim V. Then the joint distribution of Oki and Ofe^ is the same as the joint 
distribution of 

d d 
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which are therefore themselves independent. Thus we have found two hnear statistics of the r.v.'s 
Ofe that are independent. By the Skitovich-Darmois Theorem 4.2 therefore the joint distribution of 
the Ofe's is Gaussian complex valued as soon as we are able to prove that there exists at least one 
element g € G such that the real linear transformations 

C 9 z ^ Dkij{g)z and C 9 z ^ Dk2j{g)z, j — 1, . . . ,d 

are non degenerate. This follows from analyticity properties of the coefficients, as explained at the 
end of the proof of Proposition 4.8 of [2] (here connectedness is required). 

■ 
Remark that the result above does not hold for one-dimensional G-modules. As shown in [2] Example 
3.7, it is possible to construct a non-Gaussian invariant random field on the torus having all its 
coefficients independent. 

Theorem 4.1 is not really satisfactory because its assumptions are not satisfied in the case of real 
random field's, for which the coefficients are necessarily constrained by the fact that the imaginary 
parts must cancel and therefore cannot, in general, be independent. Theorem 4.1 has however its 
own interest because it contains the essence of the arguments that we use in the sequel and because 
it holds without any assumption concerning the orthonormal basis {vk)k of V . 

We consider now the case of a real G-invariant random field T. 

The results are different according to the fact that the irreducible G-module V under consideration 
is of type J^~^ or ^°, as classified at the end of section 2. 

In the case V G J^+ the fact that we deal with a real random field does not impose constraints on 
the coefficients of Ty and Ty but, of course, in order to obtain a real random field one must impose 
that in the sum Ty -\-Ty the imaginary parts cancel. In this situation the random fields will be both 
complex, in general, but their sum will give rise to a real random field. 

In this situation Theorem 4.1 states that, if dimV > 1 and the coefficients {ak)k are independent, 
then Tv and Ty are both Gaussian and their sum is a Gaussian real random field 

Conversely if F G J^'^ the coefficients must satisfy some constraints in order to ensure that the 
random field is real. It is natural then to consider the setting of an orthonormal basis that is 
self-conjugated. 

It is actually appealing to consider a basis that is formed by real functions. For such a basis, 
(wfc)fc, under the assumption that the coefficients Uk = /g2 T(x)vk{x) dx are independent. Theorem 
4.1 applies so that if the random field T is invariant then the a^'s are jointly Gaussian. Such a 
statement is not however satisfactory, as we explain below in Remark 4.6. 

We shall therefore now consider the case of a basis {vk)-i<k<i of V such that 

v^k=W (4.1) 

This means that we assume that, if the dimension of V is odd, the basis contains only one element 
Vq which is a real function. If the dimension of V is even we shall still write {vk)-t<k<t in order 
to simplify the notations (there is however no wq function). In the following arguments we shall 
consider the case where dim V is odd, the case dim V even being quite similar. 

For a basis satisfying (4.1) the fact that T is real imposes to the coefficients the requirement 
a-k = Ok- This is the usual setting in the case ^ — S"^, where, to be precise, usually one considers 
the basis of the spherical harmonics for which it holds v-k — {—l)''vk so that the condition above 
becomes a„fe = (— l)'^afe, a slight difference that does not change things. 

The argument in this case can be implemented along the same lines as in Theorem 4.1. Let us 
assume that the coefficients {ak)k>a, sue independent. Then, if mi > 0, m2 > 0, mi ^ TO2, and we 
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denote as above by Dm.m'ig) the matrix elements of the aetion of G on V , the two complex r.v.'s 

I i 

Omi = ^ Dmi,„i{g~'^)am and 0,^2 = ^ Dm2,mig~^)am (4.2) 

have the same joint distribution as ami ^^id am2 s-nd are therefore independent. The Skitovich- 
Darmois theorem cannot be applied as before, as the r.v.'s am and a_m = Wn are certainly not 
independent. But (4.2) can be written 

e 
ami = Dmi,o{g^^)aa + 2J iDmi,rn{g^^)am + Drm-m{g^^)a^^ 

rn— 1 

I 

0'm2 — Dm2,o{g )O-0 + 2_^ {Dm2,m{g )o,m + Dm2,-m{g )«« 



r?i— 1 

SO that amnO,m2 ai'G (real) linear functions of the independent r.v.'s ao, ...,«£. Therefore we can 
again apply Theorem 4.2 as soon as we prove that g (z G can be chosen so that the real linear 
applications 

Z ^ Dm,,m{g^^)z + Dmi-m{g^^)'Z, 171 ^ 1, ...,£, i ^ 1,2 (4.3) 

are non singular. It is immediate that this is equivalent to 

\Dm,M9'^)\ ^ \Dm,.-rn{g-^)\, Ul = 1, . . . , i, I = 1,2 (4.4) 

We are therefore led to state our main result under the assumption that (4.4) is fulfilled. 

Assumption 4.3 (The mixing condition) Let V C LF'{^) a self-conjugated irreducible G -module 
and let {vi)-i<i<i a self-conjugated orthonormal basis ofV let us denote by D{g) the representative 
matrix of the action of G on V . We say that (fi)-£<i<£ is mixing if there exist g (z G and < mi < 
TO2 < i such that 

l^m.,m(.g)l ^ \Dm,.-m{g)\ (4.5) 

for every < m < i, i — 1,2. 

We have therefore proved the following. 

Theorem 4.4 Assume G to be connected. Let V C L'^{5(^) be an irreducible G-module such that 
V = V. Let {vk)-e<k<e be a self-conjugated mixing (see Assumption 4-3) basis ofV. Consider the 
real random field 

Tv{x) = ^ akVk{x) (4.6) 

k 

where the r.v.'s ak,k > are independen. Then if Ty is G-invariant the r.v.'s {ak)k o,re jointly 
Gaussian and therefore also Ty is Gaussian. 

Remark 4.5 It is relevant to point out that in Theorems 4.1 and 4.4 we do not assume independence 
of the real and imaginary parts of the coefficients. Actually under this additional assumption the 
statement would become almost trivial (and much weaker) in many situations, as often the invariance 
of the random field implies that the coefficients (some of them at least, see Remark 3.5) have 
a distribution that is invariant with respect to rotations of the complex plane. And it is well 
known that this assumption together with independence of the components implies a joint Gaussian 
distribution, with no need of Assumption 4.3 (immediate consequence of the Bernstein-Kac theorem 
as recalled in Proposition 6.3 below). 
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This point is important with respect to one of the practical consequences of these results, which is 
the simulation of invariant random fields. Actually a natural and computationally efficient procedure 
in order to simulate a random field on 2^ is by sampling its Fourier coefficients. For the case 
^ = §^, for instance, Theorem 4.4 together with the fact that the basis of the spherical harmonics 
is mixing ([9] p. 145) entails that if the coefficients a^m's, to > 0, of the corresponding development 
are independent, then, in order to obtain an invariant random field, they must be Gaussian and the 
resulting random field will be Gaussian itself. Different choices of their distribution will lead to a 
random field which cannot be invariant. In particular the choice of independent r.v.'s a^m's, to. > 1 
with a complex Cauchy distribution, for example, cannot produce an invariant random field, even 
if the real and imaginary parts of aim are not independent. 

Theorem 6.2 below will imply that if the coefficients agm^s, to, > 0, of the development are 
independent then, if the resulting random field is invariant, the a^m's are Gaussian, even if the 
development is made with respect to a generic self-conjugated basis (not necessarily the spherical 
harmonics). 

Remark 4.6 If the G-module V is self-conjugated as in Theorem 4.4, it is natural to consider an 
orthonormal basis on V that is formed by real functions. For instance in the case <^ = §^ one might 
choose the orthonormal basis given by vio = Y^q and 

Vim = ^{Ytm + {-irYe,^rn), 1^^,-™ = ^ (F,™ - (-l)'"y,,_„), m>0. 

The functions u^^ are real and, if we denote a^m the coefficients of the real random field T with 
respect to the basis of the spherical harmonics, then the coefficients with respect to the basis {vim)m 
would be biQ = am and 

bim = V2Reai„i, 6^ -m = \/2Ima£^„j, m > 1, 

They are of course real r.v.'s. A repetition of the arguments of Theorem 4.1 now gives immedi- 
ately that invariance of the random field and independence of the coefhcients {h^rn)m imply joint 
Gaussianity of the coefficients {bern)m without bothering with Assumption 4.3. This would however 
be a much weaker result, as independence of the {birn)mS would imply independence of the real 
and imaginary parts of the (afm)m's, which is not required in Theorem 4.4 as pointed out above in 
Remark 4.5. 

Let V an irreducible self-conjugated G-module, Ty a real random field as in (4.6) and {vk)-e<k<e a 
self-conjugated basis as above. Then by Theorem 4.4, under Assumption 4.3, if the coefficients Uk, 
k > 0, with respect to the given basis are independent they are Gaussian. Moreover, by Corollary 
3.4, as they must have the same variance and be orthogonal, there exists c > such that for k ^ 
E[(3fiafc)^] = E[(3afe)^] — | (this is a consequence of the orthogonality of ak and a_fe — lik) and 
E[ao] = c (if the basis contains a real function vq). Conversely, is a real random field Ty with these 
properties invariant? This is the object of the next statement. 

Theorem 4.7 Let V C i^(<^) a self-conjugated irreducible G-module and {vk)-i<k<l o, self- conjugated 
orthonormal basis ofV (possibly k ^ if dim H is even). Let T a real a.s. square integrable random 
field on ^ and let (ak)-t<k<i be its random coefficients with respect to the basis above. Then if 
the real and imaginary parts of the r.v.'s ak,k > (resp. fc > if dim V is even) are centered, 
independent and Gaussian and, for k ^ 0, there exists c > such that E[(3fiafe)^] — E[(3afe)^] = | 

and E[ao] = c, then the random field 

i 



is invariant. 



k=-l 
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Proof. Wc make the proof under the assumption that dim V is odd, the case of an even dimension 
being quite similar. Therefore in the basis {vk)-i<k<i we have V-k — Vk and uq is a real function. 
Let A be the matrix of the transformation C^^+^ -^ C^^"*"^ 



Z\ 

Zq 

Z-1 

\z-,) 



/^izi + z^e) \ 



1 

V2 

1 

iV2 



(Zl 



Z-1 



zo 



(zi - Z_l) 



(4.7) 



\77l(^^"^-^V 



Lemma 4.8 below proves that the matrices D{g) = AD(g)A^^ are real orthogonal. Let a^ = Xk+iYk, 
k > 0, and ao — Z. The real r.v.'s Z, Xk, Yk, k ~ 1, . . . ,i, are independent and the matrix A maps 
the vector a = {ai, . . . , O-f)* into 

/'V2Xi\ 

V2Xi 

z 

V2Yi 

\V2YeJ 

As the distribution of a is Gaussian with all its coordinates centered and independent with a common 
variance, a is invariant in distribution under the action of every orthogonal matrix and therefore 
under the action of D{g) for every g £ G. Therefore the random vector a is invariant in distribution 
under the action of D{g) for every g £ G, which implies the invariance of Ty. 



Lemma 4.8 D{g) — AD{g)A ^ is a real orthogonal matrix for every g E G. 

Proof. It is immediate that^the rows of A are pairwise orthogonal and unitary. Therefore A is a 
unitary matrix as well as D{g). Let us prove that D{g) maps R^^+^ into M^^+^, which will end 
the proof. Let S = (^£, . . . ,^i. Cm, . . . ,77^)* G M^^+i and zu = £,k + ir/k, k > 1, z^ =^ (. Let 
z — (z£, . . . , zi, 2:0,27, • • • ,'ziY', therefore z = A^^'E. is of the form 



4^iCi+im) 



V2 



C 



h (a - ^m) 



V2 



(4.8) 



\M^^~m)l 



Now 



fe>0 *:>0 
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is a real function. Remark that the matrix A changes any vector of the form (4.8) into a vector of 
]^2£+i_ ^g D{g)z is the vector of the coefficients of the function Lgf , which is stih a real function, 
its coefficients are again of the form (4.8), so that AD{g)A^^'E. e M^^+^. 



5 On the validity of the main assumption 

In this section we investigate the validity of Assumption 4.3. 

Let us remark first that, for a given self-conjugated G-module V of L'^{^), this assumption, as 
far as we know, might be true for some orthonormal bases of V and not for other ones. So far it is 
known that it is true for the basis formed by the spherical harmonics when ^ — S'^ (see [9] p. 144, for 
a proof), if dim]/ > 3. Actually, as explained below, if dimV^ < 3 Assumption 4.3 cannot be true. 
However we investigate the implication between independence of the coefficients and Gaussianity for 
the 3-dimensional irreducible G-module of L^(§^) in Theorem 6.4. 

Remark that, as Dmk{g) ~ {LgVk,Vm), condition (4.5) is equivalent to 

\{LgVrmVmi)\ 7^ K-^g^'-m, ^'mj | for somc Q <E G and every 0<m<i?, i = l,2. (5-1) 

The main result of this section is Proposition 5.3 where we state a condition equivalent to Assumption 
4.3 carrying a more geometric meaning. This will be the key tool in the next section, where we prove 
that every self-conjugated orthonormal basis of an irreducible G-module of iy^(S^) with dim(F) > 3 
is mixing. In §7 we check the validity of Assumption 4.3 for the sphere ^ — 'S? under the action of 
G = 50(4), at least for a class of self-conjugated orthonormal bases. 
Let us first state some remarks. 

Remark 5.1 a) Mixing (Assumption 4.3), is a property of the orthonormal self-conjugated basis 
of the irreducible G-module V C i^(^) and, as far as we know, might hold for some orthonormal 
self-conjugated basis and not for other ones. However if it holds for a self-conjugated orthonormal 
basis {vk)-i<k<i, then it also holds for every other basis {wk)-e<k<e of the form Wk = Lg^Vk for 
some go <E G. Actually if D{g) denotes the matrix of the action of G on V^ with respect to the basis 
iwk)-e<k<e, that is 



Wk{g ^x) = V" Dik{g)wi{x) 



E 



-1 



then D{g) — D{gQ ggo) so that Assumption 4.3 holds also for {wk)-e<k<e- 

b) Assumption 4.3 cannot be true if Vmi is a real function of L^(^). Actually, as the left regular 
action commutes with conjugation. 



J-Jmi, — m\9) — K-'^g'^—rm'^mi) — \-l^g'VrmVmi) — \-i^gVm,V„i-) — Umi,m[9) ■ 

Therefore \Dmi,~m{g)\ — \Dmi,m{g)\ for every g (z G. This implies that Assumption 4.3 cannot be 
satisfied ii V = V and dimF = 2 or dimF = 3. Actually in the first case there is only one rrii > 0, 
whereas, if dim V = 3, the values rrii = 0,1 are possible, but vq must be a real function so that (4.5) 
cannot be satisfied for this value of nii. The case dimV = 3 is of interest because it appears in the 
Peter- Weyl decomposition of L^(§^). 

In the next section we prove however that also for this G-module if the random field is invariant 
and the coefficients ao, ai are independent, then they are necessarily Gaussian with no need of using 
the Skitovitch-Darmois theorem and therefore in a situation in which Assumption 4.3 is not satisfied. 
This raises the question whether one might prove Theorem 4.4 using a different characterization of 
the Gaussian distribution than the one provided by Theorem 4.2, possibly avoiding the need of 
Assumption 4.3 (see §8 for a more precise discussion on open questions). 
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Let H he a irreducible unitary G-modulc of real type (recall Remark 2.5) and let J : i? — > if a 
conjugation such that J^ = 1 and let us denote ( , ) the corresponding G-invariant scalar product 
on H. Let {hk)-e<k<e be a orthonormal basis of H self-conjugated with respect to J, i.e. such that 
h-k = J{hk)- We shall say that such a basis is J-mixing if (4.5) holds, now denoting by D{g) the 
matrix of the action oi G on H with respect to this basis. Of course J-mixing coincides with mixing 
if iJ C L'^i^) with the left regular action and J being the usual conjugation Jv ~ v. 

Lemma 5.2 Let H be a irreducible unitary G-module of real (resp. quaternionic) type (recall Re- 
mark 2.5) and let J : H ^^ H a conjugation such that J^ = 1 (resp. J^ = —1), then, for every 
v,w G H , 

{Jv, Jw) — {Jv, Jw) . 

Proof. It is immediate that 

{v, w)' = {Jv, Jw) 

is also a G-invariant scalar product on H, hence, by Schur lemma, there exists a real number A > 
such that {v,w)' = X{v,w). The relation J^ — 1 (resp. J^ = —1) easily implies A = 1. 

■ 
We are going to express Assumption 4.3 in terms of the action of G on the wedge product /\ H. 
Recall that /\ H is endowed with the usual G-invariant scalar product 

{Vi AWi,V2 AW2)2 := {vi,V2){wi,W2) - {vi,W2){wi,V2) . 

and we denote g{v Aw) — gu A gw the action of G on /\ H . 
Proposition 5.3 Let g G G, then |-Dmi.m(g)| = \Drni-m{g)\ if and only if 

{giKi, A h^rnJ^Kn A /l-m)2 = . 



Proof We have 



{g{hs A h^s), hm A /i-„i)2 = {ghs A gh^s, km A /i-„i)2 = 
= {ghs,hm){gh-s,h-rn) - {ghs,h^rn){gh^s,hm) = 

= {ghs, hm){gJhs, Jhra) - {ghs, h^rn){gJhs, Jh^rn) = 

- \{9hs,h^)\^ ~ \{gh^s,h,n)\^ = \D^A9)\^ - \Dm.-s{g)\^ 



where we used the fact that {Jv, Jw) ~ {Jv, Jw) thanks to Lemma 5.2. 

■ 
Assumption 4.3 can therefore be rephrased in terms of orthogonality of the G-orbits of the vectors 



2 



2 



hm A h-m in /\ H. 

To be precise, let, for every 1 < m < £, Wm C /\^ H the subspace generated by the G-orbit of 
hm A h^m] let S the set of the pairs {i,j), ^ < i,j < i, such Wi C W^-. Let S the set of the indices 
1 < i < £ such that {i,j) E S for some 1 < j < i. 

Corollary 5.4 Assumption 4-3 holds if and only if the complement set S'^ contains at least two 
indices. In particular Assumption 4.3 is verified if £ > 2 and S is empty. 

Proof. Let {i,j), I < i,j < £ and let 

F,,, -{.geG;|A,j(5)|y^|A,-,(5)|} 
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If Fi_j 7^ then it is a dense open set of G. Assumption 4.3 holds if and only if for every 1 < mi < 
1712 < ^ we have that Fmi,m 7^ for every 1 < m < £ and Fm^,„i' 7^ for every \ < m! < t. Now it 
is sufficient to observe that, by Proposition 5.3, F^j = if and only if (i, j) G S. 

■ 
In the remainder of this section we introduce a family of orthonormal bases of a G-niodule that 
arises naturally (the spherical harmonics are of this type) and for which the investigation of the 
validity of Assumption 4.3 might be simpler. 

Let H be an irreducible G-module, T a maximal torus of G and let us go back to the setting of 
Remark 3.5 and consider the decomposition (3.5). It is possible to assemble an orthonormal basis 
of H by picking a unitary vector Uk in each of the UkS. We say that such a basis is associated to 
the torus T. 

If among the C/fe's there is only one subspace at most that is associated to a given character of 
T, then the decomposition (3.5) is unique and an associated orthonormal basis is also unique, up to 
multiplication of its elements by unitary complex numbers. In this case (i.e. if among the C/fe's there 
is only one subspace at most that is associated to a given character of T) we say that H is T-simple. 

We shall see in the next sections that all irreducible sub-G-modules of i^(§^) and L^(§'^) are T 
simple with respect to the maximal tori of G = 50(3) or G = 5*0(4) respectively. 

Let us now suppose that _ff is a real G-module and denote by J a real conjugation. If -u € t/^ and 
i G T wc have, denoting u — > gu the action of G, 

tu = Xkit)u 
for some character Xk of T, so that 



t Ju — Jtu ~ Xk(t) Ju ~ x-k(t) Ju . (5-2) 

Therefore it is easy to see that an orthonormal basis of H associated to T can be chosen in such a 
way that it is self-conjugated with respect to J. We shall denote by {hk)-i<k<i such an orthonormal 
basis associated to T, where the index k ranges among the corresponding characters of T appearing 
in the decomposition (3.5). Then it is clear that if the relation 

\{ghkjhmi)\ 7^ \{gh^k,hmi)\, for some g <E G and for every fc 7^ (5-3) 

holds for the basis {hk)-e<k<e, then it holds also for every other basis that is associated to T, as 
two such bases only differ by multiplication by a unitary complex number. 

It is clear that if, in particular, H C L^(^) with the usual conjugation J and (5.3) holds, then 
also Assumption 4.3 holds for the basis {hk)-i<k<i- 

It is immediate that if H is T-simple then it is also T-simple for any other maximal torus T. 
Actually, T and T being conjugated, if T = g~^Tg, a basis {hk)~i<k<i is associated to T if and only 
if {ghk)-i<k<i is associated to T. Thanks to Remark 5.1 a), if (5.3) is satisfied for a basis associated 
to T, then it is also satisfied by all bases associated to T. 

The following result states that if an irreducible G-module H is T-simple and satisfies (5.3), then 
the same is true for every irreducible G-module that is equivalent to H. 

Proposition 5.5 Let V C -L^(^) an irreducible G-module with V = V and H a T-simple G-module 
equivalent to V. Then also V is T-simple. Moreover if (5.3) is satisfied by the orthonormal bases 
of H associated to T, then the same is true for V and every self-conjugated basis associated to a 
maximal torus of V is mixing. 

The proof is straightforward. 
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6 The sphere §^ and related examples 

In this section we prove first that, for every irreducible G-module, G = 5*0(3), of dimension > 3 
of L^(§^), every self-conjugated orthonormal basis is mixing. This extends previous results: see 
[9] p. 144 where this is proved for the basis of the spherical harmonics. We also give a proof of 
the fact that the statement of Theorem 4.4 is true for every self-conjugated orthonormal basis of 
the irreducible d-dimensional 5'0((i)-module of L^(§'^~^). This covers in particular the case of the 
three-dimensional 5'0(3)-module of L^(§^) a situation in which we know that Assumption 4.3 is not 
satisfied (Remark 5.1 b)). 

Let us recall that in the Peter- Weyl decomposition of _L^(S''), rf > 3, all the irreducible modules 
for the action of SO{d + 1) are self-conjugated (see [5] pp. 196-197), so that, when dealing with a 
real random field, in order to apply Theorem 4.4 the validity of Assumption 4.3 must be checked. 

It is well-known that 50(3) = SU{2)/{id, —id} so that the irreducible representations of 50(3) 
are the representations of SU{2) which arc trivial on {id,— id} (see again [3] or [5]). The group 
G = 5*1/(2) acts on the modules Hi formed by the homogeneous polynomials in 2 complex variables 
zi, Z2 of degree i in the following way: if p G Hg, then, if z = (zi, 22), 

gp{zi,Z2) = p{azi - bz2, bzi + az2) = p{zg) (6.1) 

where 

9^[% ^), a,6eC,H2 + |6|2 = l (6.2) 

denotes a generic element of G = SU{2). The 5C/(2)-modules Hi are irreducible and every irreducible 
5'/7(2)-module is equivalent to Hg for some ^ = 0, 1, . . . The action of —id in these representations 
is trivial if and only if i is even, so that every irreducible representation of 50(3) is equivalent to 
Hi for some £ even. 

Lemma 6.1 Let P,Q be homogeneous polynomials of degree i > I in the two complex variables 

Zi,Z2- Let 

/dP dP\ 

WQ):=det || g . (6.3) 

\ dzi dz2 1 
Then D{P, Q) is a polynomial of degree 2^—2 which vanishes if and only if P = or if Q = X P for 
some A G C. 

We give the proof of Lemma 6.1 after the following main result. 

Theorem 6.2 Let V C L'^{G), G — 50(3), an irreducible self-conjugated G-module of dimension 
2m + 1. Then, if m > 1, every self-conjugated orthonormal basis of V is mixing. 

Proof. The proof relies on the characterization of Proposition 5.3. Let £ — 2m and Hi as above. 
Then (6.3) defines a map (P, Q) -^ L){P, Q) from Hi^ Hi ^ H21-2 which is obviously bilinear and 
antisymmetric. It is also equivariant with respect to the action of SU{2) (acting both on Hi (g) Hi 
and H2£-2)- Actually, denoting dzP = (|£, |£), 

d.{gP){z) - d..P{g'-)iz) - {d.P){g'z)g' 

and one concludes easily, as detg* = 1. L{P A Q) = D{P,Q) therefore defines a linear equivariant 
map /\ Hi —>■ H21-2, such that L{P A Q) = if and only if P A Q = (thanks to Lemma 6.1). 

Let us prove that every orthonormal basis self-conjugated with respect to some conjugation J : 
Hi — > Hi (i.e. such that J{f-r) = fr) is J-mixing. Let (/_,„,...,/„) such a J-self-conjugated 
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orthonormal basis. If it were not J-mixing, then by Proposition 5.3 there would exist r, s > and 
two mutuahy orthogonal invariant subspaces Ui, U2 of /\ Hi such that fr A /_r G Ui, fs A /_s G U2 
(recall that we assume m > 1). As /^ and /_r are orthogonal, L{fr A /_r) 7^ 0, so that L does not 
vanish on Ui and by Schur lemma Ui must contain a G-submodule equivalent to i?2^-2- By the 
same argument also U2 must contain a G-submodule equivalent to H2e~2, which is not possible, as, 
by the Clebsch-Gordan decomposition (see [9] §3.5 or (7.1) below) the representation i/2£-2 appears 
only once in Hg <S> Hg and, a fortiori, in /\ Hg. 

Therefore, by Proposition 5.3, the basis {f-m, • • • , fm) is J-mixing. 

Now let (w_mj • • • J Vjn) an orthonormal self-conjugated (in the sense of ordinary conjugation, noted 
J) basis of V. The actions of G = 50(3) on V and H2m are equivalent and therefore there exists a 
map A : V ^' H2m that intertwines the two actions, that is such that ALgV = gAv for every g E G, 
V G V. Up to multiplication by a constant we can assume that A preserves the scalar product. If 
we note Jf = AJA^^f, J defines a conjugation on H2m with respect to which /,. = Avr is a self- 
conjugated orthonormal basis. By the first part of the proof we know that such a basis is J-mixing. 
Therefore there exists g G S0{3) such that 

\{LgVr,Vs)v\ = \{9frJs)H2m\ 7^ lid fr , f-s) H2m\ = \{LgVr , V^s)v\ 

and (f-rra, • • • J Vm) is mixing itself. 

■ 
Proof oi Lemma 6.1. Assume P ^ 0. If D{P,Q) = and dzP ^ 0, then there exists a function 
A : C^ ^ C such that, for every z G C, 

9,0 = X{z)d,P . (6.4) 

Recall Eulcr formula for homogeneous functions of exponent £: 

dP dP 

OZi OZ2 

and similarly for Q, so that from (6.4) we get Q = XP. On the open set C^ \ F, where P is the set 
of zeros of P, we have A = ^, so that, on C^ \ P, 

d.Q = ^d,P. 
But from A = ^ we have also, for j = 1, 2, 



^.i-fp^-o^ui-rpr^^Uo^ 



dzj P2 V^ dz-i ^ dzj ) P^V \P dzj J "^ dz-i J 

As A is analytic on C2\r, this implies that Q — const P on a non empty open set of C^ and therefore 
Q = const P everywhere. 

■ 
We address now the question of the three-dimensional irreducible G-module of L'^{EP) to which the 
previous result does not apply as Assumption 4.3 is not satisfied (see Remark 5.1 b)). 

Actually we prove a more general statement. The key argument is the following classical charac- 
terization of the normal distribution. 

Proposition 6.3 Let X — {Xi, . . . ,Xm) a W" -valued r.v. such that 

a) the distribution of X is invariant with respect to the action of SO{m); 

b) there exist i,j, 1 < i, j < m such that Xi and Xj are independent. 
Then X is Gaussian. 
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Proof. Wc can assume for simplicity that Xi,X2 are independent. Let 



A = 



1 


1 


^ 


~^ 


1 


1 


71 


71 









Vo 



... 0^ 

... 

... 

ij 



Then A G SO{d) and, by assumption, X and AX have the same distribution. In particular {Xi, X2) 
and (4= (Xi— X2), -7^ (X1+X2)) have the same distribution, so that -7= (Xi— X2) and 4= (X1+X2) 
are independent. By the classical Bernstein-Kac characterization of Gaussian measures, Xi and X2 
are therefore Gaussian (see [4] pp. 74 and 85 for a simple proof). In order to prove joint Gaussianity 
of X, just remark that rotational invariance implies that the characteristic function of X is of the 
form 



for some function ip 



i. But, as Xi is Gaussian, by choosing 6 — {9i,0, . . . ,0) we have 



V-d^l) = 0^(0) ^e-^^'^^'^c-^ -'1^1' 



where u^ = Var(Xi), which allows to conclude. 

■ 
Recall that in the Peter- Weyl decomposition of L'^{§'^^^) the smallest irreducible G- module besides 
the constants, Vd say, has always dimension d exactly (see [5] p. 197 e.g.). 

Theorem 6.4 Let (wd , . . . , u_i, wi, . . . , w_d) ford even (resp. (ui /jj_]^), . . . , wq, • • • , "^-ird-i)) /^^ 
odd d) a self-conjugated orthonormal basis, of the SO{d)-module Vd C i^(§''^^). Let T a real 
invariant random field on L'^{W^^^) such that its coefficients {ad, . . . , ai) (resp. (ai(rf_i), . • . , uq)) 
with respect to this basis are independent. Then they are Gaussian. 

Proof. Let us make the proof for d odd, d = 2m + 1. By assumption we can write a^ = X^ + ilj,, 
a_fe = Xfe — jYj. for 1 < k < d, aQ = Z , where the r.v.'s Z , {Xi, Yi),. . . ,{Xm, Y^) are independent. 

Let a = (fl-m, . . . , oo, . . . , a^)* and denote by D(g), g G G, the matrices of the left regular action 
of G = SO{d) on the G-module Vd with respect to the given orthonormal basis. By assumption the 
random vectors a and D{g)a have the same distribution for every g G G. Let now A be the matrix 
C^ — > C* defined as in (4.7) (with m replacing by t) so that 

|^/2Xm\ 
Aa= Z 

V2yi 

\^/2Y^) 

Thanks to Lemma 4.8 the matrices D{g) — AD{g)A~'^ are orthogonal. Moreover g — > D{g) is 
an irreducible representation of G of dimension d, the representations D and D being equivalent. 
As there is only one d-dimensional irreducible representation of SO{d) (up to equivalence), g — > 
A^^D{g)A is equivalent to the natural action of 5*0(6?) on C^ and therefore 1 = det{D{g)) — det D{g) 
so that D{g) G SO{d) and the image of D is SO{d) itself, as the map g — > D{g) is injective. 
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Therefore invariance of the distributfon of a with respect to the matrices D{g) entails invariance 
of the distribution of a with respect to SO{d). As the r.v.'s Xi and Z, for instance, are independent, 
Proposition 6.3 implies that a is Gaussian, and therefore also a. 

U 
The previous theorem ensures that, if T is an invariant random field on §^ and V is an irreducible 
G-module of L^(§^) of dimension 3, independence of the coefficients ai and oq with respect to any 
self-conjugated orthonormal basis of V entails Gaussianity of Ty, even if Assumption 4.3 is not true 
for such V. 

Remark 6.5 If Vd is the d-dimensional irreducible SO {d)-niodn\c and if d > 4, then Theorem 6.4 
actually follows from Theorem 4.4 as, as a consequence of Proposition 5.3, every self-conjugated 
orthonormal basis is mixing. In fact for d > 5 the module /\ Vd is irreducible whereas, for d = 4 
/\ Vd, has two irreducible components which are the eigenspaces of the Hodge* operator. Therefore 
a non zero real vector of the form ivAv cannot be contained in either eigenspace (see [3] pp. 272-274). 

Example 6.6 {S0{3) and SU{2)) In the same line of arguments it is easy to check that, for a real 
invariant random field, independence of the coefficients entails Gaussianity in the cases J^ ^ G = 
30(3) and JT = G = SU{2). 

Actually if ^ = G = SO(3) this is partially known when considering the basis given by the 
normalized columns (or rows) of the Wigner matrices: in every isotypical submodule one of the 
columns is generated by the spherical harmonics for which it is known that, for £ > 1, Assumption 
4.3 holds so that Theorem 4.4 applies. As for the other columns, they are not self-conjugated but 
conjugated pairwise, so that one can apply Theorem 4.1. 

However Theorems 6.2 and 6.4 ensure that, even considering a different decomposition of the 
isotypical spaces, it is not possible to simulate an invariant non Gaussian random field using inde- 
pendent coefficients. 

This is true also for ^ = G = SU(2) as in the Peter- Weyl decomposition, in addition to those 
already considered for G — 5*0(3), other representations appear that are quaternionic, so that the 
corresponding isotypical modules cannot contain self-conjugated irreducible modules (recall Remark 
2.5). 

7 The sphere §^ 

In this section we prove that for every irreducible G-module, G = 5*0(4), of _L^(S'^) every basis 
adapted to a maximal torus is mixing. 

We shall need some known facts about the group 50(4) and its representations. G = 50(4) is 
isomorphic to SU{2) x SU{2)/{{id,id), {—id,— id)}. Therefore its irreducible representations are of 
the form Hg (E) Hk, Hg,Hk being the irreducible modules of SU{2) introduced at the beginning of 
§6, with the condition that the action of {—id, —id) is trivial. As these modules are formed by the 
homogeneous polynomials of degree £ and k respectively in the complex variables zi, Z2, one has 

{-id, -id){p ® q) = {-l)'^+*'p®q 

and therefore the irreducible modules of 50(4) are of the form Hf, ® Hk with t~\-k even. In order to 
determine the Peter- Weyl decomposition of L^(§^), S^ = 50(4)/50(3), one must recall that, in the 
isomorphism G ~ 5/7(2) x SU{2)/{{id,id), {—id,— id)}, 50(3) is mapped into the diagonal. Then 
the action of 50(3) on Hi ® Hk is g{p<E) q) = gp<E) gq of SU{2). By the Clebsch-Gordan formula for 
SU{2), the action of SU{2) on the tensor product Hg (E) Hk can be decomposed as 

dq 

Hi(g)Hk^(^ Hi+k-2j, dq = mm{£, k) . (7.1) 
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{Ps,Pr) ^ J, -Ss.r ^ -TjrSs^r (7.2) 



Therefore the trivial representation appears in this decomposition if and only if 

. , • ^ + k ^ „ i + k ^ , 

£ + k IS even, < i, < k 

2 - ' 2 - 

that is if and only ii £ — k. We have therefore found that the representations of 50(4) appearing in 
the Peter- Weyl decomposition of L^(§'^) are exactly those that are equivalent to Hg (g) Hi. Remark 
that the smallest dimension of these, besides the case ^ = 1 of the constants, is 4, so that we do not 
have to bother with the problem of dimension 3 appearing for the sphere §^ , as discussed in Remark 
5.1 b). 

On the SU {2)-modu\e Hi introduced in §5 let us consider the polynomials Ps{zitZ2) — zlz2~''^, 
s — 0, . . . ,i which form an orthogonal basis with respect to the scalar product 

sl{ i~sy. _ J_ 

which turns out to be 5L'^(2)-invariant. Therefore the polynomials e^ ~ CgPs, s — 0, ...,i? with 

Cs = \/ [J form an orthonormal basis of the unitary SU{2)-modu\e H^. 
A maximal torus of SU{2) is the subgroup of the elements 

e-*" 
whose action on the polynomials Ps{zi, 22) = zfz2~'' is 

i,p,=e^(2-^)V. . 

Thus, with respect to the invariant scalar product (7.2), the elements Cs — CgPs with Cs — \/ () form 
an orthonormal basis of Hi that is adapted to the maximal torus T. In particular Hi is T-simple. 
The following computation is our key argument. We have 



ges = c,{azi - bz2y{hzi+az2Y-' ^ CsJ^ E ( J ( \ ) «"^"'"''fc' (-&)'"'' ^[^2""" 



r=0 h+k=r 
a<h<s 
0<k<l-s 



£ m ^ 

and therefore 

J J h+k=j \ / \ / 

0<h<s 

0<k<i-s 

Taking into account the condition h + k — j this can be written 



C.s T^ Cc 



{ges,ej) = —Hj^s ^ 

o<j-h<e-s 



-'I 

c, c, 

0<h<s 






s — h 



"^J a<h<s ^'^^ ^^ "'^ 

a<j-h<i-s 
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and therefore 

\{ges,e.j)\^ = 



1 ^ o<h<s V / V-' / / 

o<j-h<e-s ('"J) 

^'o ^ 0<h<s V / V / / 

0<j-h<t-s 

which is a homogeneous polynomial of degree 2i in the variables \a\, \b\. Let us point out that in the 
sum inside the square defining ,'3^1 the range of h is 

max(0, —1+ s + j) < h < min(s,j) . (7-4) 

Let us assume first i even, £ — 2m. Let 

Jk — ^ni+k • 

{fk)-m<k<m, IS also an orthonormal basis with respect to the 5'C/(2)-invariant scalar product (7.2) 
and adapt~ed to T. The maximal torus of SU{2) x SU{2) is T x T. As 

{te,,te,)fk^ ^ Ik, =0^(2^1^1+2^262^ /,^ ^ /,^, „^ < A;i,fc2 < m 

Xfei,fc2(*ei>*«2) 

Hi® Hi is simple with respect to Tx T and the basis {fk^ ® fk2)ki.k2 is adapted to the maximal torus 
above. Moreover if /^^ ® fk, is the eigenvector of the character Xki,k2 of T x T, then f-ki (X" f-k, 
is an eigenvector of Xj,^ j.^. We proceed now to check condition (5.3) in view of taking advantage of 
Proposition 5.5. We must show that for some mi > 0, TO2 > 

l((5i,52)(/mi «'/m2),/ri <» fr,)? 7^ l((3i,52)(/mi '^fm.,),!-,-^ (^ f-r,)]'^ for cvery - m < ri,r2<m 

(7.5) 
for some 31,52 G SU{2). We have 

|((5l,32)(/mi <»/m2),/ri <»/r2)P = l(.9l./mi , /n ) H (.g2/m2 , /r^) P 

and taking into account (7.3) 

l(5/fe,/r)P- \{9ek+m,er+m)\^ = Pk+m.r+miW\Ab\) 

SO that, denoting by ai, &i and a2, &2 the coordinates of gi and 52 in the representation (6.1), 

l((5i,52)(/™, ® /™2),/r, ® /.2)l' = ^™+™,,™+.,(|ai|, I&i|)^™"™2,™+.2(|a2|, I62I) (7.6) 

and 

\{igU92){fm,®frn,)J-r,®.f-r.)\' = ^^" ™, ,„-., (kl | , |6l |)^^"t„,,„_,, (|a2 U62 | ) • (7.7) 

In order to conclude we must prove that for some values of ai, 61, 02, 62 with |ai|2 + |5i|2 = |a2p + 
I62P — 1 the right-hand sides in (7.6) and (7.7) are different. For every r 7^ if the two polynomials 
^^ and ^m+m- m+n both homogeneous of degree Am, coincide on the circle |a|2 + |6|2 = 1, 



m+mi,m—r '^-^^'-^ =^ m-\-mi,m-\-r 
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they would coincide on the whole of M^. In order to see that this cannot happen we look at the 
monomial that exhibits the highest exponent in |a| and sec that the degrees are different. Recalling 
(7.4), we must show that the two values 

hi — min(?7i + mi,m + r) and ft,2 ~ min(m + mi,m — r) 

are different. This is done by checking directly all possibilities: as mi > 0, then 





hi 


h2 


< r < TUi 


m + TOi 


m — r 


rrii < r < m 


TO + r 


m — r 


—mi < r < 


m + r 


m + TTii 


— m < r < —mi 


m — r 


m + r 



Therefore, unless r = of course, hi ^ ft,2 in all possible occurrences. Therefore the two polynomials 
at the right-hand side of (7.6) and (7.7) arc different if one at least between ri and r2 is different 
from 0. 

Along the same lines goes the proof for £ odd. Thanks to Proposition 5.5, we have 

Theorem 7.1 Let Vi a irreducible G-module appearing in the Peter-Weyl decomposition of L'^(§^) 
of dimension > 1. Then every self-conjugated basis of Vi associated to a maximal torus is mixing. 

8 Some open questions 

This paper gives some precisions about properties of the Fourier coefScients of an invariant random 
field and clarifies some important points in the direction of characterizing the random fields on an 
homogeneous space whose coefficients in their Fourier development are independent (or at least that 
can be simulated through the generation of independent r.v.'s) on the track of [1] and [2]. 

However it also points out some natural questions that remain open to conjecture. We make here 
a tentative list. 

1) In order to prove the Gaussianity of such random field's we used the Skitovitch-Darmois the- 
orem whose application in turn requires, in many cases of interest, to ascertain that Assumption 4.3 
is verified. But we have also remarked that Gaussianity remains true in situations where Assumption 
4.3 is not true (see Remark 5.1 b). So one might think of taking advantage of a characterization of 
Gaussianity different form the one that is provided by the Skitovitch-Darmois theorem, and thus be 
ridden of Assumption 4.3. 

2) It is nevertheless of interest also to investigate the validity of Assumption 4.3. Is it always 
true (at least for real G-modules of dimension > 3)? We do not know of counterexamples so far. 

3) In a less ambitious perspective, is Assumption 4.3 true for the groups SO{d) and for the 
spheres §'^~^, d > 5? For which classes of orthonormal bases? Intuition should point towards the 
positive: as these structures contain 5*0(3) and 50(4) for which the result is proved, it might be 
possible to develop a sort of induction procedure on d in order to obtain the validity of Assumption 
4.3 or at least that the Gaussian characterization holds. Remark that one possible way of attacking 
this problem is through an extension of Proposition 6.3: does its statement remain true if the 
assumption of invariance with respect to the group SO{n) is replaced by invariance with respect to 
a subgroup of SO{n) that acts irreducibly on M"? 
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